Abstract. We have examined preliminary data of the 
globular cluster, M22, from the EPIC MOS detectors on 
board XMM-Newton. We have detected 27 X-ray sources 
within the centre of the field of view, 24 of which are new 
detections. Three sources were found within the core of 
the cluster. From spectral analysis of the X-ray sources, 
it is possible that the object at the centre of the core is a 
quiescent X-ray transient and those lying further out are 
maybe cataclysmic variables. 
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1. Introduction 

The nearby (2.6 ± 0.3 kpc, Peterson & Cudworth 1994) 
globular cluster, M22 (NGC 6656) has previously been 
studied by the X-ray satellites Rosat and Einstein (see e.g. 
Verbunt 2001 and references therein). Eight X-ray sources 
(h x < 10 34 - 5 erg s~\ Hertz & Grindlay 1983) have already 
been detected in the direction of the cluster (Johnston 
et al 1994), using the Rosat PSPC, where the source de- 
tected within the core radius is believed to be related to 
the cluster. The eight sources detected by ROSAT include 
two of the four X-ray sources detected by Einstein (Hertz 
& Grindlay 1983). In this paper we present observations 
of M22 taken with the sensitive EPIC MOS camera on 
board XMM-Newton. We have detected 27 new sources, 
using this new instrument, in the direction of the cluster. 
Through statistical and spectral analysis, we have tried to 
determine the nature of the sources detected within the 
core radius and their relationship to the cluster. 

Globular clusters are known to contain two 
classes of X-ray sources: the bright sources which 
are neutron star low-mass X-ray binaries and the 
so-called 'dim' sources, with L x < 10 34 5 erg s 1 
(Hertz & Grindlay 1983), whose nature is still dis- 
cussed. The dim sources are intrinsically faint and badly 
affected by interstellar absorption, which prevents many 
sources from being discovered. Although the nature of 
these sources is unclear, they may be cataclysmic vari- 
ables (Hertz & Grindlay 1983) or other types of binary 
systems (Verbunt & Johnston 2000). It is expected that 
globular clusters should contain many binary systems due 
to the interactions occurring within the clusters (Di Ste- 
fano & Rappaport 1994). However only a few such binary 
systems have been confirmed (e.g. Verbunt & Hasinger 
1998 and references therein). 

2. Observations and data reduction 

We have obtained approximately 37 kiloseconds (ks) of 
observations of the globular cluster M22, with XMM- 
Newton. Observations were made between 19-20 Scptem- 
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ber 2000, during the 'Routine Observing Phase'. However, 
approximately 15 ks were lost due to high background ac- 
tivity from a solar flare. We present the data obtained 
with the EPIC MOS detectors, using the full frame mode 
(see Turner et al, 2001). 

The data were processed using the standard MOS 
pipeline software, provided in Version 5.0 of the XMM- 
Newton SAS (Science Analysis Software). Bad pixels were 
removed from the data, using the task 'badpixfind' and 
the event lists were filtered, so that 0-12 of the predefined 
patterns were retained. The high background periods were 
identified by defining a count rate threshold above the low 
background rate and the periods of higher background 
counts were then eliminated from the event list. The re- 
sulting background is then very stable (see also Briel et 
al, 2001), especially on a small scale, and could then be 
subtracted. The event lists from the two MOS cameras 
were merged, to increase the signal-to-noise. The result- 
ing exposure corrected image, in the energy range 0.2-10 
keV, is presented in Figure [lj 

We compare the XMM-Newton image with the images 
from the Rosat PSPC and HRI observations, obtained 
from the Rosat archive at the Max-Planck-Institut fur ex- 
traterrestrische Physik. The PSPC data, were taken on 14 
March 1991 (8.295 ks) and the two HRI observations were 
taken on 16-17 September 1992 and 18-21 March 1993 
(9.474 ks and 31.709 ks respectively). 

3. The sources 

The XMM-Newton sources were first detected using the 
SAS EPIC source detection task 'eboxdetect', which cm- 
ploys a 'local' source detection algorithm. A box of 5 x 5 
pixels was used to detect the point sources and then the 
same box was used on the background. To detect the ex- 
tended sources, two iterations were made with a box of 10 
x 10 pixels and then 20 x 20 pixels. A detection likelihood 
could then be calculated. 

27 sources (see Table [l]) have been detected using the 
SAS task 'emldetect', using a maximum detection likeli- 
hood of 15. However, only sources in the centre of the 
field of view (FOV) have been detected with this task. 
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Fig. 1. The reduced, 22 ks exposure of M22, in the energy 
range 0.2-10keV, from the two EPIC MOS cameras. The 
inner circle (solid line) shows the core radius (85.11"). 
The larger circle (dashed line) , shows the half mass radius 
(3.3'). The right ascension and the declination are shown 
on the abscissa and the ordinate respectively. The numbers 
correspond to the ID numbers given in Table § A zoom 
of the central region can be seen in Figure 



Fig. 2. A zoom of the reduced, 22 ks exposure of M22 
(Figure |]). The inner circle (solid line) again shows the 
core radius. The larger circle (dashed line), shows the half 
mass radius. The right ascension and the declination are 
shown on the abscissa and the ordinate respectively. The 
numbers correspond to the ID numbers given in Table § 
The contours describe 3, 6, 11, 20 and 25 sigma certainties. 



We have detected a further 7 sources manually, sources 
a-g (see Table |l|), at larger off-axis angles, including the 
extended source, source c (source A in Hertz & Grindlay 
1983). These sources all have high count rates and so we 
are confident that they are real sources. All the detected 
sources can be seen, numbered, in Figure |l|. In the outer 
regions of the FOV, some background fluctuations with 
multiple contours can be seen, as fluctuations are magni- 
fied in these outer regions, when correcting for vignetting. 
It is for this reason that we have used a maximum detec- 
tion likelihood method to detect the sources. 

We have detected three sources within the core radius 
of M22 (85.11", Djorgovski 1993), sources 13, 16 and 18 
(indicated by a * in Table |l|), and five sources within the 
half-mass radius (3.3', Harris 1996), sources 13, 16, 18, 11 
and 21, where previous observations have detected only 
one (source 16: using Einstein, Hertz & Grindlay 1983; 
Rosat PSPC, Johnston et al 1994; and Rosat HRI, Verbunt 
2001). The average counts from the MOS 1 and MOS 2 
detectors with their associated errors and the unabsorbed 
flux at Earth (x 10 -14 ergs cm -2 s _1 ) in the energy band 
0.5-2.5 keV are given in Table [J This energy band has 
been chosen, so that the results are comparable with for- 
mer source detections, e.g. Johnston et al (1994). We have 
assumed a 5 keV exponential spectrum, using the column 



density Nh = 2.2 x 10 21 cm~ 2 , and a distance of 3 kpc, in 
the same way as Johnston et al (1994). The approximate 
position of each source detected and any former detection 
is also given. 

The dataset received was incomplete, so the astrome- 
try was improved using the 3 sources that were detected 
both in our data and by the Rosat HRI (sources 16, 27 and 
e, Table Q). An IDL routine was used to derive the adjust- 
ment required for the MOS data to align with the HRI 
data. We found that we required both a small transversal 
shift and rotation, in the same manner as Hasinger et al 
(2001). This leads to a residual error of ~5.5", where the 
largest error is due to the position error of the sources 
detected by the HRI. 

The six brightest sources detected by the Rosat PSPC 
have also been detected by the EPIC MOS detectors. How- 
ever, the two weakest sources detected by the PSPC are 
not apparent in our data. We have searched to a limit- 
ing luminosity of 3.0 x 10 30 ergs s _1 in the 0.5-2.5 keV 
band, derived in a similar manner as Johnston & Verbunt 
(1996), using the upper limits of 3 positions, in or near 
to the core. As the Rosat PSPC has a lower limiting lu- 
minosity of 1.3 x 10 31 ergs s^ 1 in the 0.5-2.5 keV band 
(Johnston & Verbunt 1996), but the maximum likelihood 
of existence for each is >10 (Johnston et al, 1994), the 
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Table 1. X-ray sources in the direction of M22, as deter- 
mined from the EPIC MOS observations 



For. iL4(2000) Dec(2000) Average Flux 



ID 


ID 


h m s 


O 




II 


Counts 


xl0~ 14 


1 




18 


36 


24.0 


-24 


02 


21.0 


25 ±6 


1.04 


2 




18 


36 


27.9 


-23 


44 


58.4 


47 ±7 


1.95 


3 




18 


36 


15.4 


-23 


46 


06.7 


26 ±6 


1.08 


4 


1 


18 


36 


10.5 


-23 


46 


19.6 


52 ±8 


2.16 


5 




18 


36 


34.8 


-23 


49 


19.2 


28 ±6 


1.16 


6 




18 


35 


56.8 


-23 


49 


20.6 


24 ±6 


1.00 


7 




18 


36 


40.0 


-23 


50 


48.9 


23 ±5 


0.95 


8 




18 


36 


05.9 


-23 


50 


43.1 


27 ±6 


1.12 


9 




18 


36 


15.8 


-23 


50 


51.1 


29 ±6 


1.20 


10 




18 


36 


45.7 


-23 


51 


48.0 


29 ±6 


1.20 


11 




18 


36 


18.6 


-23 


52 


09.7 


29 ±6 


1.16 


12 




18 


36 


02.7 


-23 


53 


21.2 


23 ±5 


0.95 


13* 




18 


36 


21.9 


-23 


53 


30.7 


40 ±7 


1.66 


14 




18 


36 


39.3 


—23 


53 


40.9 


43 ± 7 


1.78 


15 




18 


36 


38.7 


-23 


54 


39.5 


26 ± 6 


1.08 


16* 


R 4 


18 


36 


25.0 


-23 


54 


32.3 


193 ± 13 


8.01 


17 




18 


35 


59.1 


-23 


54 


58.4 


45 ± 7 


1.87 


18* 




18 


36 


25.2 


-23 


55 


08.8 


92 ± 9 


3.82 


19 




18 


36 


18.5 


-24 


00 


52.8 


42 ± 7 


1.74 


20 




18 


36 


03.0 


-23 


55 


21.4 


55 ± 7 


2.28 


21 




18 


36 


28.6 


-23 


56 


22.8 


23 ± 6 


0.95 


22 




18 


35 


55.6 


-23 


57 


14.1 


20 ±5 


0.83 


23 




18 


37 


00.8 


-23 


58 


24.4 


22 ±5 


0.91 


24 




18 


36 


08.5 


-23 


59 


13.9 


21 ±5 


0.87 


25 




18 


37 


00.4 


-24 


00 


22.9 


56 ±7 


2.32 


26 




18 


36 


19.3 


-23 


57 


10.0 


26 ±6 


1.08 


27 


9 


18 


36 


49.0 


-23 


48 


02.7 


42 ±9 


1.74 


a 




18 36 50.5 


-23 41 24.5 


47 ±9 


1.95 


b 




18 36 14.1 


-23 41 38.6 


64 ± 10 


2.66 


c 


A, 3 


18 35 50.6 


-23 46 49.4 


256 ± 17 


10.62 


d 




18 35 43.0 


-23 


48 


18.3 


145 ± 13 


6.02 


e 


5 


18 35 23.8 


-23 55 52.0 


171 ±16 


7.10 


/ 


7 


18 36 00.5 


-24 04 41.1 


44 ±8 


1.83 


<) 


8 


18 36 27.3 


-24 06 


48.2 


72 ±9 


3.00 



ID: 1-27, sources in the centre of the EPIC MOS FOV. 
a-g at the edge of the EPIC MOS FOV. 

For. ID: A-B, Einstein, Hertz & Grindlay 1983 
1-8, Rosat PSPC, Johnston et al 1994 
9, Rosat HRI only, Verbunt 2001 

* denotes sources within the core 



two faintest sources detected by the PSPC are likely to be 
variable sources. 

4. Cluster Membership 

It is unclear how many of the sources that we have de- 
tected are indeed members of the cluster. Searching the 
SIMBAD database, none of the X-ray sources have an op- 
tical counterpart that lies within the error box of the X-ray 



determined coordinates. We have therefore tried to eval- 
uate whether we see an overdensity of sources within the 
core of the cluster, in a similar manner to Fox et al (1996). 
We have used the log N-log S relation derived from XMM- 
Newton observations of the Lockman Hole (Hasinger et 
al 2001). Our lowest detected source has an unabsorbed 
flux value of 8.3 xl0~ 15 ergs cm" 2 s" 1 (0.5-2.5 keV). In 
the same energy range (0.5-2.0 keV) as the log N-log S 
relationship, this is 6.0 xlO -15 ergs cm -2 s _1 . Using a 
spectrum with a photon index of 2.0, the photon index 
given for sources in the energy band 0.5-2.0 keV (Hasinger 
et al 2001), rather than the 5 keV exponential spectrum 
that we have used to determine the flux, the flux drops to 
3.3 xlO -15 ergs cm~ 2 s -1 . This indicates that we should 
see 220 sources per square degree, using the log N-log S 
relationship. However the area of the core is 1.76 xlO -3 
square degrees, thus we expect 0.4 sources within the core. 
According to Poisson statistics, using the expected num- 
ber of sources and one trial, the probability of a chance 
detection of one source in the core is ^27%. 

In addition, it is also possible to calculate the prob- 
ability that the sources within the core are not simply 
spurious identifications with the cluster, in the same man- 
ner as Verbunt (2001). Using the probability, p, that one 
serendipitous source in the MOS observation is at a dis- 
tance r < R from the cluster centre, located at RA=18 /l 
36 m 24.2 s , dec=-23° 54' 12" for M22 (Djorgovski & Mey- 
lan 1993), where p = (R/rd) 2 and r<± is the radius of the 
field of view. For our observations with XMM-Newton, we 
consider only the central field of view (radius=0.134°), 
where 27 sources have been automatically detected (see 
Section |^). The probability of not finding any sources in 
a single trial simply by chance, within the core radius, 
is 97%. In 27 trials, for the 27 sources, the probability 
of finding no sources within the core radius is 44%. It is 
therefore possible that the faintest source within the core 
radius (source 13) may not be related to the cluster, where 
from their brightness, the other two core sources are more 
likely to be related. 

We also consider the sources within the half-mass ra- 
dius, as these sources are likely to be associated with the 
cluster. Both the Einstein observations, L m i„=1.02 x 10 32 
ergs s _1 in the 0.5-2.5 keV band (source D), with a lim- 
iting luminosity of the observations of approximately 1.0 
x 10 32 ergs s" 1 (Hertz & Grindlay 1983) and the Rosat 
PSPC data, where the limiting luminosity is 1.3 x 10 31 
ergs s _1 in the 0.5-2.5 keV band (Johnston et al 1994) have 
detected only 1 source within the half-mass radius. Scal- 
ing the lowest observed luminosity in a similar way as we 
have already done, and using the same log N-log S relation 
derived from XMM-Newton observations of the Lockman 
Hole (Hasinger et al 2001), for the half-mass radius of the 
Einstein observations, we should see approximately zero 
X-ray sources in such an area and similarly for the Rosat 
observations. For our observations with XMM-Newton, 2 
sources would be expected, thus we see an overdensity of 3 
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X-ray sources within the half-mass radius. If the sources 
detected within the half-mass radius are at the dis- 
tance of the cluster, their luminosities range from 
1.0-8.7x 10 31 ergs s \ in the 0.5-2.5 keV band, 
where 8.7x 10 31 ergs s 1 is the luminosity of the 
most central source, source 16. 



5. Variability 

X-ray variability can occur on different timescales, de- 
pending on the source of the variability, variations from 
milliseconds (e.g. from the neutron star spin and orbital 
motion close around the neutron star, van der Klis 2000) 
to hours (e.g. the spin period of the white dwarf in OY 
Car (Ramsay et al 2001). It can therefore be a useful di- 
agnostic of the nature of the object observed. 

We have attempted a timing study of source 16, which 
has the highest count rate of all the objects in the centre 
of the field of view. However, using the counts from all the 
energy bands, over the 37 ks, results in only 0.011 counts 
per second. We have binned the counts into bins of several 
hundred seconds and performed a period search using the 
time-series analysis package 'period'. However, we found 
no significant period in the data, other than the period on 
which the data were binned, and aliasses thereof. 

Binning the longest continuous sample of data into 
4 data bins, where each bin contains approximately 85 
counts, no variability can be seen, larger than the size of 
the error bars. We therefore require further counts to carry 
out a proper study. 



6. Spectra 

We have fitted the spectra of the sources within the cluster 
core, between 0.2-10 keV. Two of these sources are two of 
the brightest of the 27 sources detected in the centre of 
the field of view. For all of the fits, we have used a column 
density of Njj — 2.16 x 10 21 cm -2 , as was used in Johnston 
et al (1994). 

The best fit absorbed power law spectrum to the most 
central source, source 16 (Table Q) is shown in Figure |. As 
the channel boundaries are different for MOS1 and MOS2, 
the fit is duplicated, as can be seen in Figure & although 
the fit is indeed the same in both cases. The data has been 
binned up, to contain at least 20 counts in each bin. The 
fit with the lowest xt f° r the above data is a power law 
with a spectral index of 1.84±0.17 (x 2 „=0.983, 27 degrees 
of freedom). All the errors quoted are at the 90% confi- 
dence limit. However we obtain almost an equally good fit 
for a single component, bremsstrahlung model, with the 
parameters kT=4.51±1.51 keV (^=1.004). The poorest, 
yet still acceptable fit is a blackbody, the source has a tem- 
perature of 0.58±0.04 keV (^=1.406). We calculated the 
hardness ratios for the regions 3.0-10.0/1.0-3.0 keV and 
1.0-3.0/0.5-1.0 keV. These show a hard spectrum with no 




Energy (keV) 



Fig. 3. Spectrum of source 16, between 1-7 keV. Data 
from both the EPIC MOS cameras are plotted, where each 
data bin contains at least 20 counts. The solid lines show 
the best absorbed power law fit. 



high energy cutoff, with values of 0.52±0.09 and 3.05±0.99 
for the two ratios respectively. 

Sources expected are soft X-ray transients in quies- 
cence (Verbunt & Hut 1987), cataclysmic variables, RS 
CVn binaries, and recycled neutron stars (Verbunt & 
Johnston 2000). Bahcall & Lightman (1976) (see also 
Lightman et al 1980; Grindlay et al 1984) discuss the fact 
that an X-ray source in a globular cluster can be treated 
as a test particle in thermal equilibrium with the other 
cluster stars. Thus, the probability distribution for the 
location of the X-ray source depends on its mass, with 
the more massive objects lying nearer to the centre of the 
cluster. 

The most massive object expected in a globular clus- 
ter, is a quiescent X-ray transient, containing a black hole. 
In globular clusters with a high density of stars, the dy- 
namical evolution is so rapid that many black holes will 
be ejected from the system on relatively short time-scales 
(Portegies Zwart & McMillan 2000; Kulkarni et al 1993), 
although some of these will be retained in the cluster and 
may capture a 'normal star' to form a low mass X-ray 
binary, with occasional X-ray outbursts (Kulkarni et al 
1993). 

This spectrum is similar to the black hole X-ray tran- 
sient spectra depicted in Asai et al (1996), Asai et al (1998) 
and Campana & Stella (2000). If the most central object 
were a black hole X-ray transient, it would be likely that 
the compact object accretes via an advection dominated 
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accretion flow (ADAF), see Narayan et al (1996), where 
it is proposed that such objects would show hard X-ray 
spectra. 

Neutron stars have a mass of approximately 1.4Mq 
(e.g. Thorsett et al, 1993), and thus an X-ray binary in 
such a globular cluster would have a lower mass if the 
black hole were replaced with a neutron star. The ob- 
served spectrum has a blackbody temperature higher than 
that of the proposed identification of a quiescent transient 
neutron star in the globular cluster NGC 5139, and other 
similar objects examined (Rutledge et al 2001), although 
the hard spectral tail has a power law photon index that 
is consistent with quiescent transient neutron stars known 
(Rutledge et al 2001 and Asai et al 1996). Alternatively, 
globular clusters are already known to be a rich source 
of millisecond pulsars (D'Amico et al 2001). 12 millisec- 
ond pulsars (MSP) have also recently been detected by 
the X-ray telescope Chandra in the globular cluster 47 
Tuc (Grindlay et al 2000). Our spectral analysis does not 
exclude the possibility that source 16 is a MSP. 

X-ray transients have a higher minimum mass than 
cataclysmic variables, where the minimum mass of a 
known cataclysmic variable (CV) is approximately 0.5M Q 
(Ritter & Kolb 1998). A large number of CVs are expected 
to form in globular clusters through tidal capture (Hertz 
& Wood 1985) and therefore to exist and contribute to the 
number of low-luminosity X-ray sources observed in such 
clusters (Di Stefano & Rappaport 1994). RS CVn binaries 
in a globular cluster contain a late-type, usually evolved, 
(spectral type ~K/M) star. For M22, where the turn-off 
mass is about 0.8M Q (Piotto & Zoccali 1999), RS CVn 
binaries within the cluster should have a similar mass to 
CVs in the cluster. 

We have detected a second source (source 18) at the 
centre of the globular cluster, only 36" from source 16 
(source B in Hertz & Grindlay 1983). Source 16, in the 
core of the cluster, has been detected by both Einstein 
and Rosat, but in both cases, as a single source. This 
second source is, however, too faint to have been de- 
tected by the deepest survey, by the Rosat PSPC, John- 
ston et al (1994). Fitting this source in a similar man- 
ner, we find that a power law fit with a spectral index 
of 1.36±0.26 (x^=l-357, 12 degrees of freedom) gives the 
best xl- A blackbody with a temperature of 0.76±0.09 
keV (xl — 1-875), can also be used to describe the data. 
However, we find that a bremsstrahlung model is not a 
good description of the data. We calculated the hardness 
ratios in the same way as for source 16. The hardness ra- 
tios, 0.29±0.07 and 3.36±1.38 respectively, show that the 
two sources, 16 and 18, are indeed different. 

As a large number of cataclysmic variables are ex- 
pected in globular cluster systems, it could be expected 
that this source is a cataclysmic variable. However, the 
spectrum is not typical of a CV, e.g. Richman 1996, who 
finds, in general, that a brcmmstrahlung model with a 
temperature of kT=0.1-5 keV is a good description of the 



X-ray spectra of many CVs. However, Richman also finds, 
in the case of case of GQ Mus, that a low temperature 
blackbody is the best description. Van Teesling & Vcr- 
bunt (1994) find that a wider range of models fit the CV 
spectra, including multiple component fits. 

The spectra of sources 16 and 18, are significantly dif- 
ferent, especially the hardness ratios, supporting the fact 
that they are indeed different, individual sources. The 
third source in the cluster core has too few counts to be 
fitted accurately. 

7. Conclusion 

We have examined preliminary data from the EPIC MOS 
detectors on board XMM-Newton of the globular cluster 
M22. We have detected 27 X-ray sources within the field 
of view, in the line of sight to M22 and a further 7 sources 
at larger off-axis angles. These indentifications include 7 of 
the sources formerly detected, but we find no evidence for 
the two faintest sources found by the Rosat PSPC camera. 
We have detected three sources within the core of M22, 
where previous observations have detected only one. 

It is unclear which sources are associated with the clus- 
ter, as none of the sources coincide with known optically 
identified sources, although it is likely that all 3 of the 
sources detected within the core radius are associated with 
the cluster. We have presented evidence, through the X- 
ray spectra, to indicate that the X-ray source at the centre 
of the core may be a quiescent X-ray transient and that 
the two sources lying further out may be cataclysmic vari- 
ables. 

Acknowledgements. We are grateful to M. Auriere and J. 
Ballet for their comments during the preparation of this 
manuscript. We are also very grateful to the referee, Frank 
Verbunt, for many helpful comments. 

References 

Asai K., Dotani T., Mitsuda K., Hoshi R., Vaughan B., 
Tanaka Y., Inoue H., 1996, PASJ, 48, 257 

Asai K., Dotani T., Hoshi R., Tanaka Y., Robinson C.R., Ter- 
ada K., 1998, PASJ, 50, 611 

Bahcall J.N., Wolf R.A., 1976, ApJ, 209, 214 

Briel U.G., Aschenbach B., Englhauser J. et al, 1997, Rosat 
User's Handbook 

Briel U.G., Henry J.R, Lumb D.H. et al, 2001, A&A, 365, L60 

Campana S., Stella L., 2000, ApJ, 541, 849 

DAmico N., Lyne A., Manchester R.N., Possenti A., 
Camilo F., 2001, ApJ, 548, L171 

Di Stefano R., Rappaport S., 1994, ApJ, 423, 274 

Djorgovski S., Meylan G., 1993, ASP conf. series, 50, 325 

Djorgovski S., 1993, ASP conf. series, 50, 373 

Fox D., Lewin W., Margon B., van Paradijs J., Verbunt F., 
1996, MNRAS, 282, 1027 

Gray D.F., 1992, The Observation and Analysis of Stellar Pho- 
tospheres, 2nd ed. Cambridge Univ. Press, Cambridge 

Grindlay J., Hertz P., Steiner J.E., Murray S.S., Light- 
man A.P., 1984, ApJ, 282, L13 



Natalie A. Webb et al.: First XMM-Newton Observations of the Globular Cluster M22 



7 



Grindlay J., Heinke CO., Edmonds P.D., Murray S.S., 

Camilo F., 2000, A AS, 197, 126 
Harris W.E., 1996, AJ, 112, 1487 

Hasinger G., Turner J., George I., Boese G., 1992, GSFC OGIP 
Calibration Memo CAL/ROS/92-001, Technical Report, 
GSFC 

Hasinger G. et al., 2000, A&A, 365, L45 
Hertz P., Grindlay J., 1983, ApJ, 275, 105 
Hertz P., Wood K.S., 1985, ApJ, 290, 171 
Johnston H.M., Verbunt F., Hasinger G., 1994, A&A, 289, 736 
Johnston H.M., Verbunt F., 1996, A&A, 312, 80 
Kulkarni S.R., Hut P., McMillan S., 1993, Nature, 364, 421 
Lightman A. P., Hertz P., Grindlay J., 1980, ApJ, 241, 367 
Narayan R., McClintock J.E., Yi I., 1996, ApJ, 457, 821 
Peterson R.C., Cudworth K.M., 1994, ApJ, 420, 612 
Piotto G., Zoccali M., 1999, A&A, 345, 485 
Portegies Zwart S.F., McMillan S.L.W., 2000, ApJ, 528, 17 
Ramsay G., Poole T., Mason K., Cordova F., Priedhorsky W., 
Breeveld A., Much R., Osborne J., Pandel D., Potter S., 
West J., Wheatley P., 2001, A&A, 365, 288 
Richman H.R., 1996, ApJ, 462, 404 

Rutledge R.F,. Bildsten L., Brown E.F., Pavlov G., Za- 

vlin V.F., 2001, submitted to ApJ 
Ritter H., Kolb U., 1998, A&AS, 129, 83 

Thorsett S.E., Arzoumanian Z., McKinnon M.M., Taylor J.H., 

1993, ApJ, 405, 29 
Turner M.J.L., Abbey A., Arnaud M. et al, 2001, A&A, 365, 

L27 

van der Klis M., 2000, ARA&A, 38, 717 
Verbunt F., 2001, A&A, 368, 137 

Verbunt F., Bunk W., Hasinger G., Johnston H.M., 1995, 

A&A, 300, 732 
Verbunt F., Hasinger G., 1998, A&A, 336, 895 
Verbunt F., Hut P., 1987, IAUS, 125, 187 
Verbunt F., Johnston H.M., 2000, A&A, 358, 910 



